Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



A High Resolution Spectroscopic Observation of CAL 83 with 

XMM-Newton/RGS * 



o 
o 
o 

(N 
> 

o 



> 
00 

m 
o 



o 
o 



o 



X 



Frits Paerels^, Andrew P. Rasmussen^, H. W. Hartmann^, J. Heise^, A. C. Brinkman^, C. P. de Vries^, 

and J.-W. den Herder^ 

^ Columbia Astrophysics Laboratory, Columbia University, 550 West 120th St., New York, NY 10027, USA 
^ SRON Laboratory for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, the Netherlands 

Received 2 October 2000/ Accepted 27 October 2000 

Abstract. We present the first high resolution photospheric X-ray spectrum of a Supersoft X-ray Source, 
the famous CAL 83 in the Large Magellanic Cloud. The spectrum was obtained with the Reflection Grating 
Spectrometer on XMM-Newton during the Calibration/Performance Veriflcation phase of the observatory. The 
spectrum covers the range 20-40 A at an approximately constant resolution of 0.05 A, and shows very significant, 
intricate detail, that is very sensitive to the physical properties of the object. We present the results of an initial 
investigation of the spectrum, from which we draw the conclusion that the spectral structure is probably domi- 
nated by numerous absorption features due to transitions in the L-shells of the mid-Z elements and the M-shell 
of Fe, in addition to a few strong K-shell features due to CNO. 

Key words, stars: atmospheres - stars: individual (CAL 83) - white dwarfs - X-rays: stars 



1. Introduction 



Supersoft X-ray Sources (or Supersoft Sources, SSS) were 
among the first new discoveries made with ROSAT — in 
the first PSPC images of the LMC, the SSS stood out 
immediately (Triimper et al. 1991). In retrospect, some of 
these objects had been seen before with Einstein, but the 
ROSAT detection in short order of several SSS focused 
attention on them as a class. 

It is now generally accepted that the SSS are white 
dwarf stars undergoing steady nuclear burning in their 
envelopes. The inferred size of the radiating surface area, 
coupled with the characteristically very soft optically thick 
spectra, point to this idea as the natural interpretation 
for the observed emission. Some SSS are central stars 
of Planetary Nebulae, or related objects: shell-burning, 
(mostly) single stars on their way to the white dwarf 
graveyard. The classical SSS, however, is a high-luminosity 
(L ^ 0.1 iEdd)? soft {kT ^ 50 eV) low- mass binary with 
a white dwarf primary. In contrast to the case for neutron 
stars, for white dwarfs the efficiency of nuclear burning 
exceeds the efficiency of gravitational energy conversion 
by accretion, by a factor of order 100, so that SSS can 
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sustain a high luminosity with a sub-Eddington accretion 
rate. 

The basic model for the classical SSS has been worked 
out by van den Heuvel et al. (1992). A relatively high- 
mass white dwarf (0.7-1.2 Mq) can sustain steady nu- 
clear burning of accreted hydrogen near its surface if the 
accretion rate is in the range 1 — 4 x IQ^'' Mq yr^^. Such 
an accretion rate can arise if the mass donor companion 
star is more massive than the white dwarf (1.5-2.0 Mq). 
Mass transfer to the less massive star will cause the bi- 
nary orbit to shrink, and the donor will transfer mass on 
a rapid thermal timescale. The predicted luminosities are 
in the range 6 x 10'^^ — 1 x 10"^* erg s""^ , roughly in the ob- 
served range for SSS. The binary evolution aspects of this 
model, including an estimate for the birthrate and preva- 
lence of SSS binaries, have been addressed by Rappaport 
et al. (1994). 

An important part of the interest in the SSS stems 
from the fact that they may be a significant chan- 
nel to SN la events. Mass is deposited onto the white 
dwarf at such a high rate that the accreted matter can 
burn non-eruptively, which makes the white dwarf mass 
grow (in contrast to low accretion-rate binaries, such as 
classical novae, which expel most of the previously ac- 
creted matter in explosive events) , eventually reaching the 
Chandrasekhar mass. Rappaport et al. (1994) estimate 
that the rate of Galactic SN la supernovae associated with 
SSS binary evolution could be as high as ~ 0.006 yr~^. 
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Fig. 1. Exposure corrected, background subtracted spectrum of CAL 83. We show the spectrum measured with RGS 
1 (black histogram) overlaid on the spectrum measured with RGS 2 (red histogram), to emphasize the correspondence 
of coherent features in the spectrum in both spectrometer data sets. The wavelength bin equals 5 original RGS 81 ^m 
CCD spatial bins, which corresponds to approximately 0.06 A, slightly larger than one resolution element. The range 
20.1 - 23.9 A is not present in the data from RGS2, because of malfunction of the drive electronics for one CCD chip 
in this spectrometer. Included for reference are wavelengths of the strongest resonance transitions in hydrogen- and 
helium- like C. 



A high resolution soft X-ray spectroscopic observation 
can in principle directly and definitively settle any re- 
maining doubts about the nature of the emission and the 
compact object, by detecting characteristic photosphcric 
structure in the spectrum. With a high resolution spec- 
trum, we can attempt to obtain an accurate estimate for 
the fundamental stellar and binary parameters (luminos- 
ity, effective temperature, radius and mass of the white 
dwarf, mass transfer rate), which has a direct bearing on 
the evolutionary status and ultimate fate of the binary. In 
this Letter, we present preliminary results of a high resolu- 
tion spectroscopic observation of the important SSS CAL 
83 (Cowley et al. 1984, PakuU et al. 1985, Crampton et al. 
1987, Smale et al. 1988, Parmar et al. 1998) which directly 
addresses the above issues. 



2. Data Analysis 

CAL 83 was observed by XMM-Newton (Jansen et al. 
2000) during the Calibration/Performance Verification 
phase of the observatory, on April 23, 2000, for a total 
of 45.1 ksec. Data was obtained with both the Reflection 
Grating Spectrometers (RGS, den Herder et al. 2000) 
and the EPIC focal plane imaging cameras (Turner et al. 
2000). For the purpose of this paper, the high resolution 
spectrum obtained with RGS holds the primary interest, 
and we will only discuss the RGS data. 

The data were processed using custom software orig- 
inally developed for the analysis of RGS ground calibra- 
tion data, which is nearly identical in function to the RGS 
branch of the Science Analysis System (SAS). Telemetered 
CCD events are read in frame by frame and are offset- 
corrected on a pixel by pixel basis using median read- 
out maps, compiled from about 40 diagnostic images (den 
Herder et al. 2000) per CCD chip, those available from the 
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ten most recent revolutions (20 days). This process nearly 

eliminates flickering pixels from the datasct. Remaining 
flickering and warm pixels are removed from the datasets 
by thresholding duty cycle maps (for each CCD), also 
generated from about 10 contemporaneous revolutions. 
Gain and CTI correction are performed to align the sig- 
nal/energy scale across all CCD readouts. Then, event 
reconstruction is performed on a framc-by-frame basis 
where connected pixels containing significant signal are re- 
combined into composite events and the composite event 
signal is calculated from the sum of individual (corrected) 
signals. 

The standard event grade combination comprises 

events which fall within a 2 x 2 pixel region, where two 
pixels diagonally opposed to one another within the 2x2 
region are considered two separate events. The event coor- 
dinates are mapped into focal plane two-dimensional an- 
gular coordinates (dispersion and cross-dispersion). The 
dispersion coordinates are based on the known geometry of 
the RGS combined with a preliminary in-flight calibration 
based on the emission line spectra of HR 1099 and Capella 
which self-consistently determined the instrumental bore- 
sight relative to the star tracker. The event positions are 
corrected for small drifts in aspect. 

Events arc then windowed in the dispersion- 
pulscheight and focal planes, using optimum extrac- 
tion masks utilized by the response matrix generator. 
Background subtraction was performed using background 
sampling regions on cither side (in the cross dispersion di- 
rection) of the source illuminated region. Finally, a spec- 
trum file was created by histogramming the corrected 
countrates in each dispersion channel. A response gen- 
erator was used to produce an observation specific re- 
sponse matrix that provided an array of nominal wave- 
length/energy values corresponding to each channel in the 
spectrum file. This procedure was performed for both the 
first and second orders for each spectrometer, but only the 
first order spectra contained significant data. 

Figure 1 shows the exposure-corrected, background 
subtracted spectrum from RGS 1 and 2. We show the 
spectrum obtained with the two separate spectrometers 
superimposed, in order to emphasize that most of the 
observed spectral detail is real, since it appears in both 
instruments. It is immediately obvious from a glance at 
Figure 1 that blackbody spectral models, or indeed any 
smooth continuum spectrum, is completely inadequate to 
characterize the spectrum of CAL 83, and, by extrapola- 
tion, SSS in general. Previous low-resolution studies with 
ROSAT (Greiner et al. 1991; see also Heise et al. 1994) 
and BeppoSAX (Parmar et al. 1998) yielded circumstan- 
tial evidence that the spectrum had to deviate from a sim- 
ple blackbody (the bolometric correction associated with 
a blackbody yields estimates for the total luminosity that 
exceed the Eddington limit by a large factor), but the 
spectral structure that must be present in any real stellar 
photospheric spectrum is seen here for the first time. 

We have indicated the wavelengths of the strongest 
resonance lines in the H- and He-like ions of C in Figure 



1. At the approximate densities and temperatures in the 

atmosphere (corresponding to kT ^ 50 eV, log g ^ 8.5), 
CNO are expected to be primarily in their H- and He-like 
charge states, and one might expect to see their absorption 
lines in the spectrum because of their high abundances. 
Indeed, at first sight, both C VI Lya (33.73 A) and the 
C V A34.97 A n = 1 — 3 resonance line appear to line up 
with a strong absorption feature in the spectrum. Both 
the C VI Ly series limit, and the C V ionization limit 
similarly coincide with a strong feature. However, there is 
no such correspondence for any of the wavelengths in N. 
Also, closer inspection reveals the presence of other deep 
absorption features (e.g. the ones at 30.5 A, 29.0 A) that 
do not correspond to any transitions in the H- and He-like 
ions of C, N, or O. 

We illustrate this finding in the top panel of Figure 2, 
which shows the average of the two spectra from RGS 1 
and 2, compared to a model atmospheres spectrum calcu- 
lated for a hydrogen atmosphere at kTes = 45 eV, log g = 
8.5, with trace amounts of CNO added at their LMC abun- 
dance levels. This, and other models to be discussed in this 
paper, have been calculated with Hubeny's atmospheres 
code TLUSTY (Hubeny 1988; Hubeny & Lanz 1995), as 
applied to the conditions in white dwarfs (Hartmann & 
Heise 1997; Hartmann et al. 1999). We added absorption 
by neutral gas, of column density A^h = 6.5 x 10^" cm~^ 
(Gansicke et al. 1998), and convolved the model with the 
response of the RGS. Clearly, this CNO dominated spec- 
trum does not have the dense absorption structure we see 
in CAL 83. 

In our next attempt at interpretation of the spectral 
structure, we calculated a small grid of models around 
fcTeff = 45 eV, log 5 = 8.5, this time with all abundant 
elements up to Fe, at their LMC abundance levels. These 
spectra all indeed bear a qualitative resemblance to the 
spectrum of CAL 83. An example (again at feTcff = 45 
eV, \ogg = 8.5) is shown in the middle panel of Figure 

2. The spectrum has been convolved with the instrument 
response. The general structure of the model appears sim- 
ilar to the spectrum of CAL 83. But most of the structure 
visible in this model is not due to isolated, strong absorp- 
tion lines, but instead is due to large numbers of closely 
spaced, narrow absorption features from transitions in the 
L shells of Ne through Ca, as well as Fe M shell transitions. 
This dense forest of transitions is actually not resolved by 
the RGS, as is illustrated by the bottom panel in Figure 
2, where we show the same model as in the middle panel, 
but before convolution with the instrument response. 

We did not attempt a formal fit of these models to the 
measured spectrum. Even though their qualitative appear- 
ance matches the shape of the spectrum, they all differ 
from the data in significant details. In a general sense, the 
fact that the discrete absorption spectrum has not been 
resolved implies that we do not have the most powerful 
spectroscopic diagnostics that would result from the pres- 
ence of a few isolated strong transitions in single ions. The 
dependence of the shape and strength of these features 
on the properties of the atmosphere (mainly density and 
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Fig. 2. to-p -panel: Average spectrum of CAL 83 (aver- 
age of RGSl and RGS2). Superimposed is a model atmo- 
spheres spectrum calculated for /cTeg = 45 eV, loggi = 8.5, 
H CNO (LMC abundances). The spectrum of CAL 83 
clearly has a richer spectrum than this H-I-CNO atmo- 
sphere. Middle panel: Same atmosphere, with the abun- 
dant elements Ne-Ca and Fc added (LMC abundances); 
the spectrum has been convolved with the RGS instru- 
ment response. This spectrum qualitatively resembles the 
data. Lower pane};, the same spcctriim, but at the hill res- 
olution of the model calculation. The structure observed 
with the RGS is evidently not resolved. 



temperature) would have given us unambigous estimates 
of the stellar parameters. Instead, the sensitivity comes 
about through a dependence of the ionization balance of 
the mid-Z elements and Fe on temperature and density, 
and this diagnostic is of coLirsc directly subject to uncer- 
tainties in the chemical composition of the photosphere. 
The interpretation may also be sensitive to the computa- 
tional treatment of the problem and the completeness and 
accuracy of the atomic structure data used to calculate the 
line opacities, as well as to uncertainties in the treatment 
of collisional broadening of the absorption lines. We will 
address these issues in a future paper. 



3. Conclusions 

We have presented the first high resolution X-ray spec- 
trum of CAL 83. With XMM- Newton RGS, we finally 
conclusively demonstrate the photospheric origin of the 
soft X-ray emission. Moreover, our spectrum qualitatively 
resembles the spectrum expected for the emission from a 
very hot white dwarf star, and therefore confirms the basic 
model for binary Supersoft X-ray Sources. 

We detect radiation from the photosphere between 20 
and ^ 40 A, at an approximately constant wavelength 
resolution of 0.05 A. The spectrum shows very significant 
absorption structure. From comparison to model atmo- 
spheres calculations at fcTeff ~ 45 eV, logg ^ 8.5, for vary- 
ing abundances, we conclude that most of this structure is 
due to the superposition of numerous, unresolved narrow 
absorption features due to transitions in the L shells of 
the mid-Z elements, the M shell of Fe, as well as in the K 
shell of C. A quantitatively reliable estimate of the stellar 
parameters has to await a detailed analysis in terms of a 
series of dedicated model atmospheres calculations. 
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